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Abstract. A recently developed technique for the system-and-bath quantization 
of open optical cavities is applied to three resonator geometries: A one dimensional 
dielectric, a Fabry-Perot resonator, and a dielectric disk. The system-and-bath 
Hamiltonian for these geometries is derived starting from Maxwell's equations and 
employed to compute the electromagnetic fields, the resonances, and the cavity gain 
factors. Exact agreement is found with standard quantization methods based on a 
modes-of-the-universe description. Our analysis provides a microscopic justification 
for the system-and-bath quantization even in the regime of spectrally overlapping 
modes. Combined with random-matrix theory our quantization method can serve as 
a starting point for a quantum theory of wave-chaotic and disordered optical media. 

PACS numbers: 03.65.Yz, 42.50.-p, 42.50.Pq 
1. Introduction 

A standard description of losses in quantum mechanical systems is the system-and- 
bath model: One considers the system of interest coupled to a large quantum system 
acting as a reservoir or bath. The reservoir is assumed to be so large that its state 
is not much affected by the coupling to the system. In contrast, the coupling to 
the reservoir introduces both losses and noise for the system dynamics. In optical 
resonators the system-and-bath model has been used for more than 40 years [H 12 El 
as a phenomenological description for the leakage of radiation into the external 
electromagnetic field. 

In two recent papers [3J |3] we reported a microscopic derivation of the system-and- 
bath model from Maxwell's equations and derived explicit expressions for the coupling 
amplitudes between system and bath modes. We demonstrated that the model not 
only applies to good resonators with spectrally well-separated modes but that it also 
provides an exact description for spectrally-overlapping or even overdamped modes. 
The field quantization in the presence of mode overlap has generated a substantial 
literature in recent years [HUH IE] • The motivation for our work came from recent studies 
of unstable optical cavities [HI EH E] and from experiments on strongly disordered 
amplifying media O E3 , so-called random lasers. The losses in such lasers are 

typically much larger than in traditional lasers as the light is not confined by mirrors 
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Figure 1. Open optical cavities studied in this paper: (a) A one-dimensional 
dielectric slab of length I with (positive) refractive index n bounded on one side 
(x = —I) by a perfectly reflecting mirror, (b) A one-dimensional cavity bounded 
by a perfectly reflecting mirror at (x = —I) and a thin semi-transparent mirror of 
width d << I at (x = 0). (c) A two-dimensional dielectric disk with refractive index 
n. 



but by the multiple chaotic scattering within the disordered medium. The amount of 
disorder determines both the laser mode amplitudes and frequencies. Therefore, the 
modes of random media depend on the statistical properties of the underlying random 
medium, and must be analyzed in a statistical fashion. Traditional methods of field- 
quantization, such as the so-called modes-of-the-universe approach [TH] , are not suited 
for a statistical description as they do not provide explicit information about the field 
inside the resonator region. In contrast, statistics naturally enters our system-and-bath 
model |H |Hj when the internal system dynamics is modeled by random-matrix theory 
[T7] . In our previous work we have focused on the general derivation of the system- 
and-bath model, the discussion of the resulting system dynamics, and the connection 
with random-matrix theory. It is the purpose of the present paper to demonstrate our 
method explicitly for a number of models frequently used for optical resonators. 

Specifically, we consider the three types of cavities shown in figure HJ (a) A one- 
dimensional dielectric slab with (positive) refractive index n bounded on one side by a 
perfectly reflecting mirror, (b) a one-dimensional cavity defined by a perfectly reflecting 
mirror on one side and a thin semi-transparent mirror on the other side, and (c) a 
two-dimensional dielectric disk with refractive index n. The dielectrics are embedded 
in empty space. The normal modes of all three systems in question as well as their 
scattering properties can be computed exactly. Likewise, exact expressions can be 
obtained for the electromagnetic fields both within the cavities and in the external 
region. This makes all three models an ideal testing ground for our system-and-bath 
description. According to that description the Hamiltonian of the complete system, 
comprising the resonator and the external region, can be represented as 



oc 
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+ h\^2Yl duJ W Am(w) a[b m {u) + V Xm {uj) a x b m {u) + H.c. . (1) 

A m J 

There is a discrete set of cavity modes with index A, and a continuous set of external 
modes labeled by the frequency u and the discrete index m, which specifies the 
asymptotic boundary conditions far from the resonator (including polarization); u m 
is the frequency threshold above which external modes with the "channel" -index m 
exist. The operators a\ and a\ are, respectively, the bosonic annihilation and creation 
operators associated with the cavity modes. They satisfy the standard commutation 
relations [a\, ay] = 0, [a\, a y ] = 5\y. Likewise, the channel operators b m {u) and b^u) 
satisfy the commutation relations [b m (u), b n {uj')\ = 0, [b m (u), b^uj')] = 8 mn S(uj — uj'). 
Operators related to different subsystems commute. Explicit expressions for the coupling 
amplitudes W\ m (uj), V\ m {uj) will be given below. 

The electric and magnetic fields are exactly represented within the system-and- 
bath approach by 

1/2 

E(r,t)=iJ2['-^) M*)«A(r) - H.c] 




+ iJ2l duJ y— I [ b m(w,t)v m (uj,r) -H.c] , (2a) 

1/2 



B{r,t) 




a A (t)(V x u x )(r) +H.c. 



u,t)(V x v m (u,r)) + H.c. 



(26) 



where c is the velocity of light in vacuum. The modes U\(r) and v m (u, r) have support 
only in the resonator and the external region, respectively. Although the expansions (|2qj) 
and (j2fejl reduce precisely to the standard expressions known from closed systems, the 
field dynamics given by Hamiltonian (0) is fundamentally different due to the coupling 
between the resonator and the external region [3 ^] . 

The microscopic formulation of the system-and-bath model must include a 
definition of the system modes u\ and the bath modes v m (u). We will specify proper 
definitions below for the three systems of interest, but note at this point that these 
definitions (and the definition of the system-and-bath Hamiltonian itself) entail some 
degree of arbitrariness. This arbitrariness arises generically for open resonators; its 
origin is that the separation of space into two regions "inside" and "outside" the 
resonator is not unique. Consider for example a cavity with an opening, i.e. a hole in 
the material walls of the cavity boundary. Any choice of a fictitious surface covering the 
hole yields its own inside/outside separation. Moreover, different boundary conditions 
may be imposed at the chosen separating surface. Nevertheless, each such surface and 
boundary condition entail eigenmodes allowing to represent the electromagnetic field 
almost everywhere. We emphasize that the expansion cannot be expected to converge 
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uniformly or even pointwise, and in particular not on an arbitrarily chosen boundary. 
Still the system-and-bath representation is exact in the L 2 sense required by quantum 
mechanics. The freedom in the choice of the separating surface and the boundary 
conditions manifests itself in the Hamiltonian (JTJ and the field expansions (|2o|) and 
(J22J: the internal frequencies U\, the coupling amplitudes W and V, and the mode 
functions u and v, all depend explicitly on these choices. However, we show below 
that all physical observables, in particular the electromagnetic fields and the scattering 
amplitudes, do not rely on these choices and only depend on the physical boundary 
conditions imposed by Maxwell's equations. 

The outline of the paper is as follows. In section 121 we collect the main results of our 
system-and-bath approach relevant for resonant cavities. In sections the method is 
applied to the three resonators introduced earlier. We compute the respective system- 
and-bath Hamiltonians, and show that the electromagnetic fields and the scattering 
properties agree exactly with results obtained by direct solutions of the quantization 
problem. Performing the computations for different sets of boundary conditions along 
the surface separating the resonator and channel region, we demonstrate that the 
physical observables are independent of that choice of boundary condition. We conclude 
in section 

2. System— and— bath quantization 

In two recent papers |H E] we developed a microscopic system-and-bath approach to 
the field quantization in open optical cavities. To keep the present paper self-contained 
we summarize the relevant results in the present section. Our quantization technique 
applies to linear dielectric media characterized by a scalar dielectric constant e(r) that 
may explicitly depend on position and that we assume real and frequency independent. 
Working in the Coulomb gauge in the absence of charges, the exact eigenmodes f m (uj, r) 
of Maxwell's equations [12J are solutions of the Helmholtz equation 



and satisfy the generalized transversality condition V ■ [e(r)f m (uj, r)] = 0. The modes 
are labeled by the continuous frequency to and a discrete index m; the latter specifies the 
asymptotic conditions far away from the resonator. For example, these conditions could 
correspond to a scattering problem with incoming and outgoing waves. Then f m (cu,r) 
represents a solution with an incoming wave only in channel m and outgoing waves in 
all scattering channels. Using notation from scattering theory, the region outside the 
resonator will be called channel region below. We follow reference ^HJ and stick to the 
usual internal product with respect to which the differential operator in equation (j3J) is 
not Hermitian. However, the substitution 




(3) 




(4) 
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transforms equation (jSJ) into the hermitian eigenvalue problem 



Lcp r 



=V x V x 



4>m(u,r 



e r 




a; 



(5) 



for the operator L. The functions cf) m {uj) provide a complete orthonormal set in the 
space of L 2 functions. It follows from equation (J3J) that the functions f m {uj) satisfy the 
orthonormalization condition 



[dre(r)f m (u,r)f n (uj',r) = 5 mn 5(uj 



(6) 



The quantization of the electromagnetic field in terms of eigenmodes of the Helmholtz 
equation is known as the modes-of-the-universe approach [TBJ EDI EE]- Within this 
approach the field Hamiltonian reduces to a sum of independent harmonic oscillators 



1 7 

H = o E / d ^ huj (4iM4»M + A m {u)Al 



10 



and the electromagnetic fields take the standard form 



E = 1 E / dLJ 



/ da; 



'huj 

A" 

2a; y 



1/2 



1/2 



A m (a;)e- lw 7m(^r-)-H.c. 



.4, 



a; e 



(V x / m (a;,r)) + H.c. 



(7) 



(8 a) 



(86) 



The fundamental disadvantage of the modes-of-the-universe approach is that it makes 
no distinction between the cavity and the external region. As a consequence, the method 
fails to single out explicit information about the field inside the cavity. 

The system-and-bath technique is based on the separation of the electromagnetic 
fields into an inside and an outside contribution. Likewise, the eigenmodes of the 
Helmholtz equation are decomposed into contributions "living" inside the resonator 
or in the channel space. At a formal level one can achieve the inside/outside separation 
by introducing the projection operators j^U E] 



Q — dr\r)(r\ and V — / dr\r)(r 



(9) 



where \r) denotes a standard position eigenket and I represents the resonator region. 
The projection operators are orthogonal, QV = VQ = 0, and complete, Q + V = 1. 
Thus an arbitrary function in Hilbert space may be decomposed into its projections 
onto the resonator and channel space 



\4) = Q\<l>)+V\<l>) = \ii) + \v). 
Acting on \<f>) with the operator L we obtain 

L\4>) = L QQ \n) + L QV \v) + L-pglfj,) + Lpp\v), 



(10) 



where Lqq and L-p-p are the projections of L onto the resonator and channel space, 
and Lq V and L VQ the coupling terms. As explained in the introduction and illustrated 
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in references 122], this decomposition is by no means unique: different choices of 
the separating surface and different boundary conditions along that surface result in 
different decompositions of L. However, all those decompositions are subject to the 
condition that L remains Hermitian, i.e. Lqq = Lq Q , L vv = L vv and Lq V = L\,q. 

The projections Lqq and L-p-p define Hermitian eigenvalue problems for the isolated 
cavity and channel region, respectively. The closed resonator modes are the 

solutions of the eigenvalue problem 



^ssI^a) = \^fj |/*a>, (12) 

and form a discrete orthonormal basis for the cavity subspace. Likewise, the channels 
modes \u m (u)) satisfy the equation 

L PV \u m {uS)) = \v m (u)), (13) 

and constitute a continuous basis for the channel region. We note that the functions 
H\{r) have support only inside the resonator, while the functions u m (u,r) are nonzero 
only in the channel region. 

It follows from these definitions that the resonator and channel modes form a 
complete set of Hilbert space functions in the respective subregions. One may therefore 
use these functions to expand the eigenstates of the Helmholtz equation 

oo 

<f> m (uj,r) = J2 a m\(uj)nx(r) + J duJ ' Pmm?(w, u')v m >(u', r), (14a) 

A m ' <*> m > 

oo 

f m (uj,r) =J2 a m\(u)u x (r) + J2 J <ko' (3 mml (uj,uj')v m r(u',r), (146) 



2 



where we introduced the functions u\(r) = ii\(r)/Je{r) and v m (u>,r) = 



v m (uj, r)jyj e(r). The expansion coefficients can be recovered from the mode functions 
using the relations 

a m \(u}) = (fj, x \4>m(u)) , (15a) 
P mm >(u),u') = (v m ,(uj')\(j> m (uj)). (156) 

The system-and-bath Hamiltonian ((JJ and the field representations ()2a)l and (J26j) 
follow upon the introduction of bosonic creation and annihilation operators a and for 
the cavity modes, and similar operators b and W for the channel modes. Finally, the 
coupling amplitudes in the Hamiltonian (JIJ are of the form 

c 2 

W\ m (uj) = - (p\\LQp\i/ m (u)) , (16a) 

c 2 

V Xm (to) = (fj.* x \Lop\u m (u)). (166) 

2 v /^7 

The notation \fi*) means (r\[i*) = /x*(r). For time reversal invariant systems the wave 
functions may be chosen real, then the amplitudes W and V become real and identical, 
W = V. 



Quantum theory of multimode fields: Applications to optical resonators 



7 



In order to test the system-and-bath approach for the resonators of interest we 
compute several physical observables below. Examples include the electromagnetic 
fields, and the eigenmodes of the Helmholtz equation. We also study the cavity 
resonances, i.e. the complex frequencies that determine the cavity response to external 
excitations in the presence of the coupling to the outside world. Formally the resonances 
are found as the poles of the resolvent operator Q(u) projected onto the cavity space and 
analytically continued in the second Riemann sheet After projection the resolvent 
can be written as 

Gqq(u) = -, x2 » ( 17 ) 

,7 -LM 



where the non-Hermitian operator L e s(u>) is expressed through the projections of the 
differential operator L, 

L c s(u) = Lq Q + Lq V 2 L VQ , 

[fj -L vv + ie 

u m {u')){u m {u')\ 



Lqq + c L Qr ^ / du L-pq, (1 

m J u — 10 + le 



here the limit e —>■ + is implied. The second line follows upon using the completeness 
of the channels modes with the help of equation (|T3|) . The operator L e ff(o;) is closely 
related to the effective Hamiltonian in open quantum systems |22j . To determine the 
system resonances one must solve the eigenvalue problem, 

LeffCuO&H) = er?(w)|&(w)). (19) 

Since L eS (u) depends parametrically on u>, both its right eigenstates iCi^)) an d the 
complex eigenvalues Oi{u)) generally depend on uj as well. The states |£i(u;)) correspond 
to the Kapur-Peierls states [2U EE| of scattering theory. Combining equations (fTTj) 
and ()19|) the cavity resonances are found as the solutions of the fixed point equation 
co~i{uS) = uj. We show below that this resonance condition is independent of the 
inside/outside separation and the choice of boundary condition made within the system- 
and-bath description. 

A well-known quantity that shows resonant behavior is the cavity gain factor 

J reI drp (u,r) 

that is closely related to the dwell time of scattered radiation inside the resonator [26j . 
Here, p(u, r) is the local density of states inside the cavity and po(u,r) the free space 
local density of states in the absence of the cavity. The cavity gain factor thus measures 
the change in the local density of states introduced by the cavity. For later use we 
note that the integrated density of states can be expressed in terms of the expansions 
coefficients (I15aj) . 

/ drp(^,r)=££K A H| 2 . (21) 
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3. One dimensional dielectric cavity 

Our first example is the one dimensional dielectric cavity depicted in figure ^a) |27j . 
The dielectric with refractive index n is nonabsorbing and nondispersive. It is bounded 
by a perfectly reflecting mirror at x = — I while there is no mirror at the other end of 
the dielectric at x = 0. The free space outside the cavity runs from x = to infinity, 
and light propagates freely there, n = 1. We assume the electromagnetic field to be 
linearly polarized with the electric field vector pointing in the ^-direction. Cavity field 
excitations will decay due to leakage into the empty half-space. The dielectric function 
of the total system including the cavity and the attached half-space is given by 

e(x) =n 2 Q(-x) +0(x), (22) 

where the Heavysi de-function Q(x) is equal to one for positive x and vanishes for 
negative x. The exact eigenmodes of Maxwell's equations for this problem are given 
in equation (jA.lj) . To solve the problem within the system-and-bath approach, we 
separate system and bath at x — 0. The cavity thus runs from x = — I to x = 0, and 
the channel region from x = to oo. The boundary conditions at the interface are only 
restricted by the requirement that they lead to an Hermitian eigenvalue problem. Below 
we address two different such boundary conditions: In the first case we set Neumann 
boundary conditions for the cavity; Hermiticity [22 El then imposes Dirichlet conditions 
for the channel problem. In the second case we consider the inverse situation with 
Dirichlet boundary condition for the cavity and Neumann conditions outside. 



3.1. Cavity with von-Neumann boundary conditions 



The inside/outside decomposition of the differential operator L reads [3] 

1 d 2 , 5(x-0„) d 
fi{x) + 



Lqq/J,(x) 

Lppv{x) 
L VQ fi(x) 
Lqpv[x) 



n 2 dx 2 
d 2 



n- 



dx 



dx 2 
S'(x - 



v{x) -5'(x-0 + )v(0+), 



n 



S(x 



^(0- 
0_) d 



-v{x) 



(23a) 

(236) 
(23 c) 
(23cf) 



The shorthands T indicate the limits where the interface at x 



n dx' 

is approached from 

inside respectively outside the resonator. The singular terms guarantee the matching 
conditions for the electromagnetic field at the interface. In addition, these terms ensure 
the Hermiticity of Lqq and L-p-p. The range of the operators in equations ()23ajl and 
(123 6|) within Hilbert space is given by the functions for which the singular term vanishes. 

The eigenmodes [i\{x) of the closed cavity are the solutions of the eigenvalue 
problem Lqq/z a (x) = k 2 fi\(x) with k\ = uj\/c. From equation (|23aj) . they satisfy 
the equation 

d 2 

— nx(x) + n 2 k 2 x ^x(x) = 0, (24) 
dx A 
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subject to the boundary conditions 



fi X (-l) = 0, (25a) 

^-f*\(x) = 0. (256) 

The first condition is imposed by the perfectly reflecting mirror at x = — I, and the 
second follows from the requirement that the singular boundary term applied to fi\ 
must vanish. The normalized solutions of the eigenvalue problem form the discrete set 

12 

fJ>x(x) = y y sm(nk\(x + /)), (26) 

with wave numbers k\ = (2A + l)7r/2ni (A = 0,1,2,...). In the channel region, the 
eigenvalue problem reads 
d 2 



Mk, x) + k 2 u(k, x) = 0, (k = oj/c), (27) 
ax 2 

with Dirichlet conditions at the resonator surface, 

u(k,0+) = 0. (28) 
This determines a continuous set of (^-normalized channel modes, 

[2 

v(k, x) = \ — sm(kx). (29) 

V 71 

Since both the cavity and channel modes are real valued functions, the coupling 
amplitudes W and V become real and identical, W = V. Combining equations (jl6aj) . 
()23a| with the mode functions (J26|) and (|29|) . we obtain 



m(k) = Vx(k)= [ -^J — . (30) 

n V Ttk\l 

The result for the internal frequencies uj\ = ck\, together with the coupling amplitudes 
(j30|) . and the mode functions (J26|) and (J29|) . completely specify the system-and-bath 
Hamiltonian and the electric and magnetic field. 

We now turn to an illustration of our results. To compare the exact scattering states 
with their representation in terms of cavity and channel modes, we combine equations 
(HHI, (flBljl with the results (TP) . (|2l)jl . and pi. This yields the expansion coefficients 

(-l) x+1 I k kcos(nkl) 

a\{k) = - n — -^7—, (31) 

nyixlyk 1 — kj) 

(3(k, k') = ±-\v ( 2k '^^ 2 k) ) - Ml - S k )6(k' - k)} . (32) 



The symbol V denotes the principal value. The S-matrix S k and the amplitude Ik 
are given, respectively, by equations (jA.2|) and (jA.3|) . Figure 121 shows the real part of 
the scattering wave function with wavenumber kl = 18. We compare the exact result 
(solid gray line) with the system-and-bath expansion (dashed line). In the resonator 
region we only included 11 cavity modes with wavenumber centered around kl = 18. 
The agreement is very good; deviations are only visible close to x — 0, i.e. near the 
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Figure 2. Real part of the scattering wave function corresponding to kl — 18, 
computed for a one dimensional dielectric cavity with refractive index n — 1.5. The 
solid line is the exact solution, the dashed line the expansion in terms of the resonator 
and channels modes. Only 11 cavity modes with k\l around kl = 18 were included. 
The inset shows the cavity gain factor as function of kl for a range around kl = 18. 



boundary separating system and bath. It has been argued before |28j that cavity 
or channel expansions must fail close to the boundary; so a remark concerning the 
status of such expansions is in order here: The inclusion of all cavity and all channel 
modes yields an exact point-to-point representation of the scattering function and its 
derivative, everywhere exception the point x = 0. This representation does not converge 
uniformly but it is exact in the L 2 sense. Therefore the system-and-bath expansion is 
an exact representation of the scattering state in the underlying Hilbert space. 

To determine the cavity resonances we solve the eigenvalue problem for the non- 
Hermitian operator L e ^{k). Explicit calculation ( Appendix B ) shows that L e g(k) acts 
on an arbitrary resonator state fi(x) according to 



L eS (k)n(x) 



1 d 2 



n- 



dx 2 



6(x-0. 



n- 



dx 



— ikp(0- 



(33) 



Due to the singular term, the action of L eS (k) generally goes beyond Hilbert space. 
The range of L e s(k) within Hilbert space is defined by the wave functions for which the 
singular term vanishes. It follows that the right eigenstates £j(k,x) are solutions of the 
Helmholtz equation 

d 2 

-^(k,x) +n 2 a 2 (k)^j(k,x) 



dx 2 



0. 



that obey the boundary conditions 
-0 = 0, 



-^•(fc, x , 



dx- 



(34) 

(35a) 
(356) 



The first condition results from the perfect mirror at x = — I while the second defines 
the so-called Siegert boundary condition. It accounts for the leakage out of the cavity: 
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In the channel region that boundary condition implies a purely outgoing wave. For a 
fixed value of k, one finds the discrete set of solutions 

£j(k,x) = Aj(k) sm(no-j(k)(x + /)), (36) 

with some normalization factors Aj(k). Substituting the solutions into equation (j35fc|) 
we obtain the secular equation for the eigenvalues, 



CTj(k) = —ktan(ncrj(k)l). 



n 



(37) 



The fixed point equation k = o~j(k) determines the cavity resonances. Analytical 
continuation of equation (J37|) into the complex plane then yields the resonance condition 

tan(nW) + in = 0, (38) 

which has solutions only for complex k. The resonances can be found analytically and 
are given by 



1 

nl 



(2j + 1)tt , i 



2 ln (M); 



3 



0,1, 



j7r + -ln(|r|); j 



1,2, 



(n> 1), 
(n< 1), 



(39) 



where r = (n— 1) / (n+1) is the reflection amplitude at the dielectric surface. Comparison 
with the direct calculation (cf. equation (jA.2)0 shows that the resonances coincide with 
the poles of the scattering-matrix. All resonances have the same width and are located 
along a straight line in the lower half of the complex plane. The resonance spacing, 
i.e. the difference in real parts of two successive resonances, is constant, A = ixjnl. 
The resonances start to overlap when the modulus of the reflection amplitude becomes 
smaller than \r\ = exp(— 7r). 

We finally evaluate the cavity gain factor. The free-space local density of states is 
Po(k,x) = sin 2 (A;(x + /)). Integration over the cavity volume yields 
o 



/ dxp (k,x) = - 

J 71 



1 - 



sin(2£;Z) 
2kl 



(40) 



The integrated cavity density of states follows from equations (|2*Tj) and (|3~Tj) by means 
of the Poisson sum rule, 



dx p(k, x) 



\Ik\ 2 k 2 cos 2 (nkl) 



txIw 



E 



i\h 



4tt 



A=0 

sin(2nkl) 
2nkl 



2\2 



(41) 



Combining these results with equation 
G c {k) 



n 2 


-I sin(2nfci) 
2nkl 




n 2 cos 2 (nkl) + sin 2 (nfc/) 




1 - 


sin(2M)" 
2kl 



2J) for Ik, we obtain the cavity gain factor 

(42) 



The inset in figure 121 shows the cavity gain factor over a range of kl. The peaks are 
equally spaced and have approximately the same high and width as expected from 
equation (jH9*|) . 
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3.2. Cavity with Dirichlet boundary conditions 

It is interesting to carry out the system-and-bath quantization in a basis other than that 
considered in the previous section. To that end we reconsider the dielectric resonator of 
figure E^a) but perform the system-and-bath quantization with interchanged boundary 
conditions at the resonator /channel interface: The resonator modes are now required to 
satisfy Dirichlet boundary conditions at x = while Neumann conditions hold for the 
channel modes. The differential operators corresponding to this choice have the form 

1 d 2 8'(x - 0_) 
L QQ fi(x) = -rf-j-jK x ) + ~ 2 M0-) , (43a) 



d 2 , x „ , d 



x'=Q- 



(436) 

=0+ 

(43c) 



6(x-0+) d , 

Lp Q fi{x) = -—fx(x 

n ax 

5'(x - ) 

L QV v{x) = — — f(0+). {43d) 
n 

The closed cavity eigenmodes of Lqq solve the Helmholtz equation and satisfy Dirichlet 
boundary conditions both at x = I and x = 0_. The second of these conditions 
follows from the requirement that the application of Lqq on any eigenmode must yield 
a vanishing singular contribution. The eigenmodes form the discrete set of functions 

H\{x) = y|sin(n/c A (x + I)), (44) 

with eigenvalues k\ = nX/nl (A = 1, 2, . . .). In a similar fashion one finds the continuous 
set of channel modes 

[2 

v(k, x) = J — cos(kx), (45) 

that satisfy Neumann boundary conditions at x = 0+. Substituting the mode functions 
into the definitions (|16aj) and (jl6fe|) we obtain the coupling amplitudes 




^,=V„ = (-l)^. (46) 

We note that the cavity eigenfrequencies and the coupling amplitudes obtained with 
the present set of boundary conditions differ from the results obtained in the previous 
section. Consequently, two different system-and-bath Hamiltonians are obtained in the 
two cases. However, as we show below both Hamiltonians provide an equivalent, and 
exact, description of the field dynamics. 

Expanding the modes-of-the-universe f(k,x) in terms of resonator and channel 
modes, we find the expansions coefficients 



ny irl(k 2 

(1 

\ k' 2 




V ] - i7r(l + S k )5(k' - k) 



(4£ 
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cavity with the same parameters as in figure The solid gray curve is the exact 
solution. The system-and-bath expansion (dotted, dashed line) is based on 11, 
respectively, 25 cavity modes satisfying Dirichlet boundary conditions at x = 0. The 
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with S*fc and denned, respectively, by equations (jA.2j) and (jA.3|) . In figure 0] we 
compare the exact scattering wave function ()A.1|) (solid gray line) with the mode 
expansion in terms of the expansion coefficients (}4Tj) and (|48)1 . Perfect agreement is 
found in the channel region x > 0. In the cavity region there is slow convergence close 
to x = 0, due to the Dirichlet boundary conditions at x = 0_. The slower convergence 
visible in figure El must be compared with the faster convergence found for the other set 
of boundary conditions (Figure EJ). It indicates that, in spite of the freedom inherent 
in the projection formalism for the choice of boundary conditions, certain boundary 
conditions are better suited for the problem yielding good approximations with less 
terms in the mode expansions. 

As in section 13.11 we can now evaluate the system resonances and the cavity gain 
factor. With the present boundary conditions the operator L e g(k) reduces to (see 
|Appendix B ) 



L eS (k)/i(x) 



1 d 2 



:fx(x) 



6'(x-0. 



n- 



k ax 



x'=0- 



.(49) 



n 2 dx 2 

It is illustrative to compare this with the result (J3"3")l that holds for interchanged 
resonator/channel boundary conditions. Both results differ in their singular terms. 
However, upon projection onto the Hilbert space the same operator is recovered as the 
singular contributions vanish. In both cases the resulting boundary condition at the 
resonator /channel-interface is the Siegert condition (J35&J) . 

The integrated cavity density of states follows upon combination of equations (|2Tj) 
and (jUJ), with the result 
o 



dx p(k, x) 



l\h\ 



1 - 



sin(2nkl) 
2nkl 



(50) 
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It agrees with the result (|41j) obtained for the other set of boundary conditions. This 
demonstrates that the physical observables are indeed independent of the choice of 
boundary conditions. 



4. One dimensional cavity with a semitransparent mirror 



The model of Ley and Loudon (29J is a one dimensional cavity defined by a totally 
reflecting mirror at one end and a semitransparent mirror at the other end (Figure ^b)). 
The electric field is linearly polarized in the z-direction. Radiation can leak out through 
the semitransparent mirror modeled by a dielectric slab of width d and refractive index 
n. The limit d — > and n — > oo with n 2 d = rj fixed is taken at the end of the calculation, 
here rj is a factor characterizing the mirror transparency. In this limit, the frequency 
dependent mirror reflection and transmission amplitudes are given by 

^ = 2^5" m = 2^' (51) 
They obey the common relations for symmetric mirrors, |r| 2 + \t\ 2 = 1 and rt* + r*t = 0. 

The exact eigenmodes of Maxwell's equations for this problem are given in 
equation (|A.4|) . Within the system-and-bath approach there are two natural ways 
of a resonator /channel separation: Either one assumes the mirror to be part of the 
cavity or the mirror is part of the channel region. Here, we stick to the latter choice. 
Accordingly, the cavity runs from x = — I to x = and the channel region from x = 
to oo. The alternative definition with the mirror being part of the cavity can easily 
be shown to lead to the same physical results. We choose Dirichlet conditions for the 
resonator boundary at x = 0_, which implies von Neumann conditions for the outside 
problem. The differential operators corresponding to these definitions are 

d 2 

Lqqh(x) = ——^n(x) + 5 (x - 0_)/i(0_), (52a) 



L-p-pv(x) 



dx 2 

1 d 2 fv(x)\ 5{x-0 



n(x) dx 2 \n(x) 




(526) 

x'=0+ 

(52 c) 



5(x-0+) d 

L ^ {x) = -^xT^ {x ' 

Lopv{x) =- < J'(x-0_)4^4, (52d) 

n(0+) 

where n(x) is the refractive index in the channel region, 

<x) = + " ~ h (53) 

[1 (d<x). 

The eigenvalue problem defined by Lqq reduces to that of the dielectric resonator of 
section EOl in the case when the dielectric function equals 1. Adopting our earlier results 
in that limiting case we find the closed cavity eigenmodes 

Hx(x) = Jjsin(k x (x + L)), (54) 
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with the eigenvalues k x 
the Helmholtz equation 

d 2 



(A = 1,2,...). The channel modes are the solutions of 



u(k,x)+n (x)k u(k,x) 



with Neumann boundary conditions at x 
conditions 



—v(k, d- 
n 

l_d_ 

n dx 



u(k, d A 



u(k, x) 



dx 



u(k, x) 



■ 0, (55) 
+ . In addition, they must satisfy the two 

(56) 
(57) 



imposed by the continuity of the electric and magnetic field at the right end of the 
semitransparent mirror. The shorthands d± indicate the limit where d is approached 
from the left (oL) or from the right (d + ). Solving for v(k, x) and taking the limit d — > 0, 
n — > oo with n 2 d = r], one obtains the following continuous set of channel modes, 

1 ' -ikx i c iv.\Jkx\ o n.\ '^ — ^k 



u(k, x) 



+ S c (k)e 



S c (k) 



(5f 



"2-k % i + rjk 

The coupling amplitudes follow upon substituting the wavefunctions (jHljl . (|?)Hj) into the 
definitions ()16aj) and ()16fe|) . The result is 



W x (k) = V x (k) 




(59) 



1 — irjk 

Finally, the representation of the exact modes f(k,x) in terms of the system and bath 
modes yields the expansion coefficients 

[—l) x Ikk\ sm(nkl) 



a x {k) = 
P(k,k r ) 



7Tl{k 2 - k\) 



1 

2^ 



V 



l-S*(k')S k , S*(k') - S k 



m(l + S*(k)S k )5(k' - k) 



(60) 



(61) 



k' — k k' + k 

where Sk and Ik are given by equation (|A.5|) and equation (jA.6|) . respectively. Figure |U 
shows the real part of the scattering wave function with kl = 28.9. The exact solution 
(IA.4j) (solid gray line) is compared with the mode expansion (dashed line) using the 
expansion coefficients ()60|) and (|61)1. The first 35 cavity modes were included. The 
deviations from the exact scattering wave function visible near x = can be made 
arbitrary small by including more terms in the mode expansion. 

In order to find the resonances one must solve the eigenvalue equation for L e s(k). 
The real space representation of L e s(k) for our choice of boundary conditions follows 
upon combination of equations (fTHj) and (j52qj) (see Appendix B| ), with the result 



L eG (k)n(x) 



dx- 



:fi(x) + S'(x - 0. 



M0-) - 



i 



(62) 



k(i + r)k) dx' 

where fi(x) is an arbitrary resonator state. The range of L e g(k) within Hilbert space is 
defined by the resonator functions for which the singular term on the right hand side 
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Figure 4. Real part of the scattering wave function with fcZ = 28.9 for a 
one dimensional optical cavity with a perfectly reflecting mirror at i = — I and 
a semitransparent mirror at x = 0. The mirror transparency is characterized by 
7/ = 0.0453 corresponding to the reflection coefficient \r(kl = 28.9) | 2 = 0.3. The solid 
line is the exact solution, the dashed line represents the mode expansion truncated to 
the first 35 modes in the cavity region. Inset: Cavity gain factor as function of kl for 
a range around kl = 28. 



vanishes. In particular, this holds for the right eigenstates £j(k, x). It follows that these 
states satisfy the boundary condition 

1 d 



x=0_ 



(63) 



k(i + rjk) dx 

at the semitransparent mirror. There is a discrete set of solutions, 

£j(k,x) = Aj(k) sin(a j (A;)(a; + L)), (64) 

with some normalization constant Aj(k). Substitution into equation (}63|) yields the 
equation for the eigenvalues o~j(k), 

o-j(k) cot(aj(k)l) = k(i + krj). (65) 

After analytical continuation of the fixed point equation k = o~j(k) into the complex 
plane we obtain the resonance condition 

i + r]k - cot(fcZ) = , (66) 

that coincides with the equation for the poles of the S*-matrix (cf. equation (|A.5J) ). 

To quantify the resonant response of the cavity to external excitations we compute 
the integrated local density of states, again using the Poisson sum rule, 

l\I k \ 2 sm 2 {kl) +2? k 2 x 



dx p(k, x) 



An 



1 - 



sm(2kl)~ 
2kl 



(67) 
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Combination with the free-space density of states (J4(J|) gives the cavity gain factor 

1 



G c (k) 



(68) 



1 — rjk sin(2/c/) + (nk) 2 sm 2 (nkl) 

With increasing kl sharper resonances are found in the cavity gain factor (see 
figure Efinset)). The reason is the reduction of the mirror transmission for large M 
that, in turn, enhances the lifetime of the cavity resonances. 



5. Dielectric disk 



In this section we demonstrate our quantization technique for resonators of spatial 
dimension larger than one. Specifically, we consider a two-dimensional circular dielectric 
of radius R and refractive index n (see figure EJc)). The resonator is embedded in free 
space. We restrict ourselves to TM modes with the electric field polarized in the z- 
direction. It is convenient to use polar coordinates r = (r, <fi) below. The dielectric 
function then reads 

e (r) = n 2 Q(R-r) + Q(r- R). (69) 

The scattering problem at the resonator can be solved exactly. The exact eigenmodes 
of Maxwell's equations are summarized in equation (jA.7|) . 

To apply our quantization technique we separate system and bath along the 
boundary of the dielectric: The dielectric disk (r < is taken as the cavity, while 
the free space (r > R + ) becomes the channel region. For the cavity we assume Dirichlet 
boundary conditions at r = which implies Neumann conditions at r = R + for the 
channel problem. The differential operator resulting from this choice reads 



1 d 

L<2QV(r, <p) = ^V 2 /i(r-, 0) + 

n 



L PV v(r, 0) 
L VQ fi(r } 0) 
L QV u(r,(f)) 



-V 2 z/(r, 
8(r - i? H 



dr 
5(r - R 



-5(r-R-)) %M^- 
r / rr 



d 



n 



dr 



7^(^,0) 



I — 

dr 



-v{r\ 



■— (-8(r 
dr \r 



R. 



—u(R, 
n 



(70 a) 
(706) 

(70c) 

(70 d) 



Due to the rotational symmetry we can choose the eigenstates to be angular momentum 
eigenstates. 

The eigenmodes fi m \ of the closed cavity are labeled by the angular momentum 
number m and the radial quantum number A. They solve the Helmholtz equation 

V 2 /i mA (r, 0) = -n 2 /4 A /i mA (r, 0), (71) 

and satisfy the Dirichlet condition /i mA (i?_,0) = 0. The normalized eigenstates are 
given in terms of Bessel functions of the first kind, 



AVaO", 0) 



<fifRJ m+1 (x m x) 



(72) 
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with m = 0, ±1, ±2, . . . and A = 0, 1, 2, . . .. The eigenvalues are k m \ = x m \/nR where 
x m \ denotes the A-th zero of J m {r). In a similar fashion, one determines the eigenstates 
in the channel region. They can be written in terms of Hankel functions, 

Vm{KrA) = Xhf^WKkr) + S m (k)H£\kr)), (73) 

V 07T 

with the diagonal element of the scattering matrix 

S m (k) = - ™ 1 J . (74) 
H$\kR) 

The channel states obey the Neumann condition -j^p m (k,r,(j)) r =0. 

In addition to the internal frequencies uj m \ = ck m \ we need the coupling amplitudes 
W and V to fully determine the system-and-bath Hamiltonian. Combination of 
equation (j70dj) with the mode functions (|T2"j) . (fT3~j) and the definitions (jl6aj) and (jl66|) 
yields after a short calculation 



W m x,n( k ) = - p 5 m,n, (75) 

TikRHm (kR) 



nkRHm (kR) 

We note that the resonant amplitude W couples only cavity and channel modes with 
the same angular momentum, while the antiresonant amplitude V couples modes with 
opposite angular momentum. This feature guarantees angular momentum conservation: 
The resonant terms Wxma^bm and W^ m a ?T) ,6j n account for the creation of a photon 
with angular momentum m and the simultaneous annihilation of a second photon with 
the same angular momentum. By contrast, the antiresonant terms, V\mCL m \b_ m and 
V\m a m\b- m describe the simultaneous annihilation or creation of two photons with 
opposite value of angular momentum. In both cases, the total angular momentum 
is conserved. 

We now turn to the electromagnetic field and the cavity resonances. In the cavity 
region the exact scattering states f m {k) can be represented in terms of the cavity modes 
u m \, with the expansion coefficients 

/,\ _ k kmxImhJm^nkR) (77\ 

a m xA«) ~ y 2 n(k*-k 2 mX ) m ' n} [ ] 

where I m k is given by equation (jA.9|) . It suffices to compare the radial component of 
the scattering wave functions. In figure we show the real part of that component for 
angular momentum m = 13 and kR = 10.5. The solid gray line is the exact result (|A.7|) 
while the dotted and dashed lines represent the system-and-bath expansion using the 
first 11 and 25 cavity modes, respectively. 

The cavity resonances are obtained by solving the eigenvalue problem for L e ^{k). 
The calculation is presented in | Appendix B| and yields the resonance condition 

J m (nkR)H^\kR) - nJ' m (nkR)H£>(kR) = 0, (78) 
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Figure 5. Real part of the radial component of the scattering wave function with 
angular momentum m = 13 and kR — 10.5 inside a dielectric disk with radius R and 
index of refraction n = 3.3. The solid gray curve is the exact result. The dotted 
(dashed) line follows from the mode expansion taking into account the first 11 (25) 
cavity modes. 



which is equivalent to the equation that determines the poles of the S'-matrix (cf. 
equation (|A.8[) ). 

Substitution of equation (|77jl into equation pTf yields the integrated local density 
of states inside the dielectric disk (see Appendix Q , 

/+oo l p2 It 12 
drp(k,r)= ]T J ^ (j 2 m (nkR) - J m+1 (nkR)J m ^(nkR)) . (79) 

disk m - 00 

Together with the free-space local density of states po(k) = k/2n, we obtain the cavity 
gain factor 

= 4:(Jj(nkR) - J m+1 (nkR)J m -i(nkR)) 

(nnkR) 2 \J m (nkR)H'£\kR) - nJ ! m (nkR)H { rk\kR)f 

where we used the explicit expression (IA.9j) for the mode amplitude I m k. The cavity 
gain factor displays a set of very sharp resonances (see figure EJ), corresponding to states 
with angular momentum kR < m < nkR, superimposed over a smooth background due 
to broad resonances with m < kR [30J. 



6. Conclusion 



In this paper we have derived exact system-and-bath Hamiltonians for a number of 
optical resonators. Physical observables computed with these Hamiltonian agree with 
the results obtained by direct scattering methods. It follows that the system-and- 
bath approach, originally developed as a phenomenological model, can provide an exact 
quantitative description of optical systems even in the regime of overlapping modes. 

On a formal level our quantization method is based on the Feshbach projector 
technique [21] that has extensively been used in nuclear and condensed matter physics 
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Figure 6. Cavity gain factor as function of kR for a dielectric disk with radius R and 
index of refraction n — 3.3. The sharp resonances correspond to states with angular 
momentum 10 < m < 18, very sharp resonances with m> 18 are not resolved. 



[HU E2j • In that areas the method has become a powerful tool for the description of 
disordered and chaotic media. Statistical theories for such media have been obtained 
employing a random-matrix assumption for the underlying system Hamiltonian |17j . 
Our application of the projector technique to the field of quantum optics may lay the 
ground for a statistical treatment of disordered and wave-chaotic optical media. 
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Appendix A. Exact modes of Maxwell's equations 

In this appendix we summarize the exact solutions of Maxwell's equations for the 
three systems treated in the paper; for a more detailed derivation we refer to references 
[2*71 [23 HJU]. The modes-of-the-universe are taken to be scattering states with an 
incoming wave in only one scattering channel. In the 2d example, this channel is labeled 
by the angular momentum number m. The scattering states are normalized according 
to equation ©. 

Appendix A.l. One dimensional dielectric cavity 

For the dielectric cavity of figure GJa) the scattering states are given by 




(-1 < x < 0), 
(0 < x). 



(A.l) 
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They satisfy the boundary condition f(k, —I) = 0, imposed by the completely reflecting 
mirror at x — — I, and are continuous with continuous derivative for any value of x > —I. 
The single-channel S-matrix and the mode strength amplitude are given, respectively, 
by 

n + itan(nfc/) 
n — itan(nfc/J 
—2m 

(A.3) 



ncos{nkl) — isin(n&;/) 
Appendix A. 2. Cavity with semitransparent mirror 

The scattering states for the one dimensional cavity with the semitransparent mirror 
(Fig. db)) have the form 

1 [ I k sia(k(x + I)) (-1 < x < 0), 

f(k jX ) — I (A.4) 
y/2n 1 exp(— ikx) + Sk exp(ikx) (0 < x), 

where the scattering matrix S^. and the mode strength amplitude Ik are given by 

Sk = i - yk + cotjkl) ^ 

i + rjk — cot(kl) ' 

2i IK s 

(A.6) 



(i + krf) sin(/c/) — cos(kl) 

Here n specifies the mirror transparency. The modes satisfy f(k, —I) = at the perfectly 
reflecting mirror and are continuous everywhere else. At the semitransparent mirror 
their derivative has a discontinuity proportional to the mode amplitude, f'(k, 0+) — 
/'(&;, 0_) = —r]k 2 f(k,0), where the prime denotes differentiation with respect to the 
position. 

Appendix A.3. Dielectric disk 

The exact eigenstates for a two dimensional dielectric disk of radius R and refractive 
index n embedded in empty space read 

Ik ■ j. I ImkJm(nkr) (0 < r < R), 

f m (k,r,(j) ) = \— e- im(t> { 1, 1 m V (A.7 

1 } VSvr \H^\kr) + S mk H^\kr) (R<r), 1 J 

where m labels angular momentum. The channel with index m is open when k exceeds 
the channel threshold k m = m/nR. Due to rotational symmetry angular momentum is 
conserved, and the S'-matrix is diagonal in the angular momentum basis 

_ H'M(kR) - n[J' m {nkR)/J m {nkR)}Hj){kR) 

fJmk — 77T\ 7T\ * I ^/ 

Hl 1} (kR) - n[J' m {nkR)/J m {nkR)}H { rk\kR) 
The mode strength amplitude takes the form 

4i 



irkR (j m {nkR)H'£\kR) - nJ' m (nkR)H$ \kR 



(A.9) 
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Appendix B. External Green functions 

Here we evaluate the action of the non-Hermitian differential operator L e g(k) on an 
arbitrary cavity state. According to equation (fTSj) L eff (fc) is the sum of two operators. 
The first contribution Lqq has already been computed in equations (|23ajl . (|43aJ) . (j52aj) 
and (j7()oj) . The second contribution has the form 

L QV g ch (k)L VQ ^r,<P) = Lqp £ Jdk' ^^^J ' Wfr 0), (B.l) 

n /-J-). 

where £/ c h stands for the retarded Green function of the isolated channel region. Below 
we compute this contribution for the cavities of interest. 



Appendix B.l. One dimensional channel 



For a one dimensional semi-infinite channel x > with Dirichlet boundary conditions 
at x = + we find the retarded Green function from the solutions (J29j) by contour 
integration, 



/ ^ f 

Q ch (k,x,x) = - dk 

7T J 



hi' 



f sin(k'x) sm(k'x') 
k 2 - k' 2 + ie 



ife|x+x'| —ik\x—x'\ 



2k 



(B.2) 



Combination with the definitions ()23cj) . (|23d|l yields the real space representation of 
L e ff(k). The relevant derivatives have to be done with care as the limits x — > and 



x 



must be taken independently. The result takes the form 



L Q vGch{k)L VQ n(x) 



iks{x - o_; 



n- 



(B.3) 



Together with equation ()23aj) we arrive at the result fpTTT^> for L e g(k). 

In a similar fashion we evaluate the retarded Green function for the isolated channel 
problem with Neumann conditions at x = 0+. Using the solutions (f43J) we obtain in 
this case 



Gch(k, x, x') 



71 



dk 



,cos(k'x) cos(fcV) 



L (Jk\x+x'\ 

?,k \ 



k 2 - k' 2 + ie 

I Q -ik\x-x'\ 



2k 

Combination with the real space representation 
respectively, yields 

t n fu\T / \ iS'(x-O-) d 



(B.4) 

and (|43(ijl of L V q and L Q -p, 



n 2 k 



dx r 



x'=0- 



(B.5) 



Combinations of equations (fTHj) . ()43aj) . and (|B.5j) yields L eS (k) given in equation (}4T?j) . 
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The semitransparent mirror has width d and refractive index n. In the limit d — > 0, 
n — > oo with n 2 d = 77 fixed, the retarded Green function can be expressed in terms of 
an integral over products of the scattering states (f5S|). 



Gch(k, x, x') = — j 



dk 



k 2 - k' 2 + ie 



(B.6) 



The integral can be done by contour integration. Using the unitarity of S c (k) and taking 
into account that S c (k) is analytic in the upper half of the complex plane, equation (jB.6|) 
reduces to 

Gch(k,x,x') = -± (e ife ^^'l + S c (k)e ik \* +X '{ 

Combination with the definitions (J52c|) and (|52dj) yields 

5'(x-0J) d 



L Q vGch(k)L PQ fi(x) 



k(i + rjk) dx 



x'=0- 



(B.7) 
(B.8) 



Equations 



B.8|) along with equation (JTHj) yield equation (J52j) . 



Appendix B.3. Angular momentum channels 

The retarded Green function for the two-dimensional Helmholtz equation with Neumann 
boundary conditions along a disk of radius R, has the form 



G ch (k,r',(j)',r,< 



i •^ c irn(<j>-(f ) ') Hm (kr~ 



Hm\kR) 

x (HW(k r< )H'W(kR) - H'i 2 \kR)H£\k r< )) , (B.9) 
where r< (r>) stands for the smaller (larger) of r and r'. Substitution of the definitions 



flZig) and (l70l yields 

LQvGch(k)L-p Q n(r,(j)) = 



R d 
2n 2 n dr 



^5(r-R-.j\ 



x 



+00 

E 



2tt 



m=— 00 



, eim( ^ 0') 







kH2\kR)dr 

Combination with equations (fTHjl and 
resonator state /i(r, (p) like 

n 

1 „ 2?r 



= R- 



(B.10) 



L eS (k)n(r,(f)) 



V 2 fJL(r,<j>) + ^-£-(-6(r-R- 
n z or \r 



shows that L e ff(k) acts on an arbitrary 

n(R-,<f>) 



1 -1-00 

— y 

2vr ^ 



m=— 00 



i' e im( 



(B.ll) 



kH'£\kR)dr'i 

Conservation of angular momentum and the requirement that the singular terms must 
vanishes for the right eigenstates £ m j(k,r,(f)) yield the boundary condition 

HW(kR) d 



Cmj [ki R- 



kH'X\kR)9r 



„tmj(k,r",(p) 



--R- 



(B.12) 
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There is a discrete set of solutions 

UAk,r,(/>) = A mj (ky m tj m (na mj (k)r) } (B.13) 

with normalization constants A m j(k). Substitution into equation (|B.12J) yields the 
equation for the eigenvalues 

kJ m (na mj (k)R)H^(kR) - na mj (k)J' m (na mj (k)R)Hg\kR) = 0. (B.14) 

Combination with the fixed point equation k = a m j(k) finally gives the resonance 
condition (|78|). 

Appendix C. Local density of states inside the dielectric disk 

Here we determine the density of states inside a dielectric disk with radius R and 
refractive index n. We start from equation (|21jl. Using the coefficients (|77jl. we obtain 

+00 7,1 r |2 00 p 

drp(k,r)= £ k ^Jl{nkR)Y, ( J m t, , 2 • (CI) 



2n 2 mv f^ikz-k 



m\) 



disk 

We concentrate on the evaluation of the last sum on the right hand side. It can be 
written derivative 



kL, Id™ 1 



E ^mA _ _____ \ ^ 



(C.2) 



The right hand side can be simplified by noticing that nRk m \ = x m \ is the A-th zero of 
the Bessel function J m (x). The sum on the right hand side can then be identify as the 
trace of the Green function for the radial part of the Helmholtz equation, in a dielectric 
disk with Dirichlet conditions on the disk perimeter. Equation (jC.2|) becomes 

00 k 2 1 B r 

E (ie%ji = -.s^rtWV^.ol^. (C.3) 

The disk Green function can be found by standard methods [33J and is given by 

YK 

Gdisk(k/a, r,r) = J rn (nkr < /a) 
A.J m (nkR/a) 

x (Jm (nkR/a) Y m (nkr > /a) — Y m (nkR/a) J m (nkr > /a)) , (C.4) 

where Y m is a Bessel function of the second kind and r< (r>) the smaller (larger) of r 
and r'. The problem then reduces to the evaluation of the integral in equation ()G3|) . 
Using the relations 

R 2 



JdrrJ 2 m (kr) = — (j 2 m (kR) - J m+1 (kR)J m ^(kR)) , (C.5) 


r R 2 1 

JdrrJ m (kr)Y m (kr) = —(j m (kR)Y m (kR) - -J m+1 (kR)Y m _ 1 (kR) 



-\ J m -i(kR)Y m+1 (kR)) , (C.6) 
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we find after some straightforward manipulations the result 

R 

/ rrg ^< fc/a - r - r)= 2^/ a) - (a7) 
Substitution into equation (|C.3j) yields 

oo k2 p2 

E (p _ m ^ A)2 = 4^P^) i J ^ nkR ) - J m+ i(nkR)J m ^(nkR)) . (C.8) 
After substitution into equation (jC.l|) one arrives at the integrated density of states 
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